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t-H ' ABSTRACT 

o : 

' Context. The question how the initial conditions in a star-forming region affect the resulting mass function of the forming stars is one 

' of the most fundamental open topics in star formation theory. 

| Aims. We want to characterize the properties of the cold dust clumps in the Carina Nebula Complex, which is one of the most massive 

Qstar forming regions in our Galaxy and shows a very high level of massive star feedback. We derive the Clump Mass Function (C1MF), 
explore the reliability of different clump extraction algorithms, and investigate the influence of the temperatures within the clouds on 
QQ ' the resulting shape of the C1MF. 

y—( | Methods. We analyze a 1.25° x 1.25° wide-field submillimeter map obtained with LABOCA at the APEX telescope, which pro- 

vides the first spatially complete survey of the clouds in the Carina Nebula Complex. We use the three clump-finding algorithms 
i i ■ CLUMPFIND, GAUSSCLUMPS and SExtractor to identify individual clumps and determine their total fluxes. In addition to as- 

' suming a common "typical" temperature for all clouds, we also employ an empirical relation between cloud column densities and 

O temperature to determine an estimate of the individual clump temperatures, and use this to determine individual clump masses. 
Results. We find that the CIMFs resulting from the different extraction methods show considerable differences in their shape. While 
(-H . the C1MF based on the CLUMPFIND extraction is very well described by a power-law (for clump masses well above the completeness 

limit), the CIMFs based on the extractions with GAUSSCLUMPS and SExtractor are better represented by a log-normal distribution. 
I 1 We also find that the use of individual clump temperatures leads to a shallower C1MF slope than the (often used) assumption of a 

O ■ common temperature (e.g. 20 K) of all clumps. 

jH ■ Conclusions. The power-law of dN/dM oc M~ 195 we find for the CLUMPFIND sample is in good agreement with C1MF slopes found 

in previous studies of the CIMFs of other regions. The dependence of the C1MF shape (power-law versus log-normal distribution) on 
the employed extraction method suggests that observational determinations of the C1MF shape yields only very limited information 
about the true structure of the cloud. Interpretations of log-normal C1MF shape as a signature of turbulent pre-stellar clouds versus 
power-law CIMFs as a signature of star-forming clouds may be taken with caution for a single extraction algorithm without additional 

^ ' information. 

ON ' Key words. Stars: formation - ISM: clouds - ISM: structure - ISM: individual objects: NGC 3372 - Submillimeter: ISM 

o\ • 
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^» ! 1. Introduction Despite these differences in the formation environment, the 

mass function of the forming stars, i.e. the final result of the 

^~ j ; The process of star formation occurs in a wide variety of envi- star formation proc ess, appears to be remarka bly unifo r m (see 

\ ronments and corresponding physical conditions. Most nearby iBastian et al] l2010l and references therein). ISalpeterl (ll955l) 

. . ■ (d £ 300 pc) star-forming regions are low-density clusters or showed that the distribution of initial stellar masses can be 

>- associations in which only low- and intermediate mass stars described by a power-law dN/dM oc M~ a with a = 2.35. 

form.The interaction between the young stars in such regions Further studies confirmed the power-law shape of the Initial 

is minimal and they can thus be considered as forming essen- Mass Function (IMF) with a power-law index a b etween 2 1 and 

^ \ tially in isolation. The more massive and generally more dis- 2.5 in the upper mass regime (above ~ 1Mb ) dKroupall20"0lt 

■ tant star-forming regions, on the other hand, contain high-mass ISchneider & Brooksl 120041: lAndre et all 1201 Ol) and showed a 

stars (M > 20 M ), too. These massive stars profoundly in- shallower slope and a turn-over for lower masses. Alte rnatively 

fluence their environments by creating HII regions, generating the IM F can be described by a log-normal distribution (IChabrierl 

wind-blown bubbles, and, finally, exploding as supernovae. This |2003l) . 

massive star feedback can disperse the natal molecular clouds, e . A , . , . • . e 

F ,. , , , , , One of the most fundamental open questions in star forma- 

but ionization fronts and expanding superbubbles can also com- tiofl ± ifj how ^ MF ^ mM tQ ^ ^ molec _ 

press nearby clouds and may thereby trigger the formation of ^ doud denfji stmctme Submillimeter observations have 

new generations of stars. , • . . . , c ^ .■ , , ■ 

& become an important tool for investigating these dense cloud 

structures. LABOCA has already been used for surveys of dis- 

tant m olecular clouds and clumps (Sc huller et al . 2009; Bo t et al.l 

* Based on data acquired with the Atacama Pathfinder Experiment hni nK 

(APEX). APEX is a collaboration between the Max-Planck-Institut fiir 

Radioastronomie, the European Southern Observatory, and the Onsala In observations of the generally distant massive star-forming 
Space Observatory. regions, the individual cloud cores can usually not be resolved at 
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submillimeter or radio wavelengths. The relevant structures ac- 
cessible to the observations are then the clumps within a cloud 
and t he corresponding Clump Mass Function (C1MF). Molecular 
line dKrameret all 119981: IWong et all 120081) and dust con tin- 
uum emission (I Johnstone & B ally 2006T Munozetal ]|2003) ob- 
servations of several star-forming regions show that the C1MF 
of molecular clouds also can be described by a power-law 
distribution. The slope of th e C1MF, with a ~ 1.4-2.0 
dElmegreen & Falgaronefl996l) . is, however, typically shallower 
than the core mass function and the stellar IMF. 

Recent studies of molecular cloud structure suggest that 
the C1MF or the distribution function of the column den- 
sity (N-PDF) can be used as an indicator for the evolution- 
ary state of a molecular cloud dKainulainen et al.ll2009l 1201 it 
iBallesteros-Paredes et al.ll201 ll) . In clouds in which the star for- 
mation process has not yet started, turbulence is expected to lead 
to a log-normal distribution. As soon as the star formation pro- 
cess starts, the denser structures get dominated by gravity and 
this should result in a power-law distribution of the masses (or 
column densities) in the upper mass (density) range. Therefore, 
the shape of the observed C1MF or density distribution function 
is sometimes considered as an indicator for the physical status 
of a cloud. One problem with such an interpretation is that dif- 
ferent observing and analysis techniques s ometimes yield differ- 
ent C1 MF shapes for the same cloud (see lReid & Wilsonll2006t 
iReid et aHl2010t) . Further determinations of the C1MF of clouds 
in star-forming regions are required to gain more insight into 
these questions. 

The Carina Nebula Complex represe nts one of the most mas- 
sive Galactic star-forming regions (see ISmith&Brooksll2008l 
for an overview). It is located in the Sagittarius-Carina spiral 
arm at a very well known distance of 2.3 kpc (Smith 2002), 
and hosts at least 65 O stars (in cluding several Q3 star s and the 
02 supergiant HD 93129A; see lSmifh & Brooksll2007l) . several 
Wolf-Rayet stars, and the Luminous Blue Variable n Carinae, 
which is our Galaxy's most luminous star known and expected 
to explode as a supernova within the next Myr. These massive 
sta rs are located in several loose clusters (Tr 14, Tr 15, Tr 16, 
see iTrumplerll 1 9 3 OT) . with ages ranging from around 1 to several 
Myr. The feedback from the massive stars has already dispersed 
most of the initial clouds in the central region. The radiation and 
stellar winds of the massive st ars have form ed numerous giant 
dust pillars (South Pillars, see ISmifhll2006l) located a few par- 
sees away from the stars in south-eastern and north-western di- 
rection. Several studies found clear indications f or ongoing and 
triggered star formation in the tips of these pillars (tMegeath et al.l 
fl996llSmifh & Brooksll2008l ISmith et al.ll2010bl) . 

The stellar populations in the Carina Nebula Complex have 
recently been studied in detail in the context of comprehensive 
multi- wavelength surveys. A deep wide-field Chandra X-ray 
surveys has allow ed to detect a large sample of the low-mass 
stellar population dTownslev et al.ll2.01 it iPreibisch et al.ll201 lat 
IWang et alj|201 It IWolk et alj|201 ll) . A very deep near-infrared 
survey of the central region with HAWK-I at the ESO Very 
Large Telescope allowed a detailed char acterization of the prop - 
erties of the X-ray selected young stars (IPreibisch et al.ll201 lbl) . 
Additional information about the recent star formation processes 
came from Hubble Space Telesc ope (HST) observations of pro- 
tostellar jets dSmifh et al.l2010aj) and Spitz.er infrared imaging of 
the South Pillars region dSmifh et"al]|2010bl) . 

While the stellar populations are by now rather well inves- 
tigated, the available observations of the molecular clouds (and 
thus on the interaction of the massive stars and the surround- 
ing clouds) were quite limited until recently. The most extensive 



existing data set of the Carina Nebula Complex at radio wave- 
lengths is a N ANTEN survey i n seve ral CO lines, covering a 
4° x 2° area by lYonekura et al.l (120051 ). Mopra 12 CO (1-0) data 
of a smaller part of the central r egion, but wi t h high er spatial 
resolution has been presen ted bv lBrooks et al.l (fl998) (see also 
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ISchneider & Brooksl 

We have recently performed a large-scale (1.25° x 1.25°) 
submillimeter mapping of the Carina Nebula Complex with 
LABOCA at the APEX telescope to obtain de tailed informa- 
tion on the structure of the cold dusty clouds (Pr eibisch et al.l 
1201 id) . Our LABOCA data revealed the very clumpy structure 
of the clouds. We found that the total mass of the dense clouds to 
which LABOCA is sensitive is ~ 60 000 M . This value agrees 
fairly well with the mass estimates for the well localized molec- 
ular gas traced by 13 CO (lYonekura et al. 2005). This high mass 
emphasizes that despite several mega-years of ongoing cloud de- 
struction due to massive star feedback, there is still a very large 
amount of cloud material available for future star formation. 

In this paper, we analyse the distribution of the cloud 
masses in the LABOCA map. We us e the three commo n 
clump-finding algor ithms CLUMPFIND dWilliams et alj|1994l) . 
GAUSSCLUMPS dStutzki & Guestenl 1 19901) and ^Extractor 
dBertin & Arnoutsl 19961) to identify individual cloud clumps and 
to measure their total sum-mm fluxes (Sect. 2). In Sect. 3 we 
estimate individual temperatures for each clump and calculate 
the clump masses. We then discuss our results and the follow- 
ing consequences. Finally, Sect. 4 summarizes the results of this 
study. 



2. Observation and data analysis 

2.1. Observations 

To investigate the clumpy dust structures and to determine the 
C1MF in the Carina Nebula we used the sub-mm emission data 
of the 1.25° x 1.25° wide-field map we derived from observa- 
tions we perf ormed at the 1 2-me ter APEX ("Atacama Pathfinder 
Experiment". iGiisten et al.l2006l) tel escope, with the Larg e Apex 
BOlometer CAmera LABOCA (see ISiringo et ai1l2009l) . which 
operates in the atmospheric window at 870 fim. LABOCA has an 
angular resolution of 18.6" which, at the distance of the Carina 
Nebula, allows us to resolve structures down to 0.2 pc. The map 
thus is suitable to resolve the molecular clumps, but not the in- 
dividual cores and provides the first spatially complete survey 
of the dust clouds in the Carina Nebula. These clumps are ex- 
pected to be the sites at which star clusters can form, contrary 
to prestellar cores which are suppo sed to form single or gravita- 
tionally bound multiple protostars dWilliams et al.ll200"ob . 

For the d ata reduction the BO lometer array Analysis soft- 
ware (BOA) dSchuller et al.l l2009h has been used, resulting in 
a pixel size of 6.07" in the final map. The surface brightness 
can be transformed to integrated fluxes by multiplying it with 
the pixel-to-beam-size ratio (0.0941 beams/pixel). The average 
rms noise level for the map is about 20 mJy/beam, which corre- 
sponds for isolated compact clumps with estimated temperatures 
of T =s 20 - 30 K to a sensitivity limit in mass of about 2M . 
For the clouds we measure intensities up to around 4 Jy/beam. 
The total flux measured in the map above a 3<x noise level 
amounts 1147 Jy. The observations are described in detail in 
IPreibisch et"alld201 Id) . 
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Fig. 1. A 20 arcminutes wide detail from our LABOCA map at 870 fim of the Tr 14 region in the Carina Nebula Complex with the 
positions of the clumps found with the different clump-finding algorithms. The units of the scale bar on the right are Jy/beam. In the 
left panel the CLUMPFIND sample is shown with circles of the effective radial size of the clumps calculated by the algorithm. The 
GAUSSCLUMPS (middle panel) and SExtractor (right panel) samples are shown as ellipses. The sizes of their axes correspond to 
the rms deviations. The bright source in the lower left corner is the luminous blue variable r\ Car. 
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Fig. 2. The flux distribution of clumps for the three different 
clump-finding algorithms with a statistical error -J jj^? for the 
CLUMPFIND sample. 

2.2. Clump-finding algorithms 
2.2.1. CLUMPFIND 

The CLUMPFIND (CF) algorithm from IWilliams et all ( fl994T) 
contours the data at a given threshold, which is recommended to 
be chosen as three times the rms noise of the map. Afterwards 
it identifies the peaks of emission which locate the clumps. It 
then contours the data in discrete steps defined by the user down 
to lower intensities until the lowest threshold is reached, assign- 
ing the pixels within the contours to the before located maxima. 
Blended structures are assigned via a friends-to-friends algo- 
rithm to the closest peak. No specific clump profile is assumed. 

As we find an rms noise level of about 20 mJy/beam for our 
data we start contouring at an intensity level of 60 mJy/beam and 
go in steps of 40 mJy/beam up to 4.5 Jy/beam, below which all 
the emission of the clouds is found. CLUMPFIND extracts 687 
molecular clumps with a total flux of 720 Jy from the LABOCA 
map. This is about 63% of the total flux in the map. A section of 
the central region around the star cluster Tr 14, with the clumps 



detected by CLUMPFIND marked as blue circles, is shown in 
the left panel of Fig. [1] The strong point-like source 77 Car has 
been excluded from this and the following samples. Fig.|2]shows 
the distribution of the extracted total fluxes F tot . 

2.2.2. GAUSSCLUMPS 

GAUSSCLUMPS (GC) dStutzki & Guested 1 19901) is a least 
square fitting algorithm, that assumes the clumps to have a 
Gaussian shape and works directly on the continuous intensity 
distribution of the data. It first searches for the highest intensity 
peak on the data map and fits a Gaussian shaped clump around 
it. Then the algorithm subtracts the clump from the original map 
generating a new map on which it continues with the next itera- 
tion step. 

The GAUSSCLUMPS algorithm calculates an rms noise 
level of 22 mJy/beam, and we define a detection threshold of 
3cr, so we can compare the results. 371 clumps are detected, 
only about half the number found by CLUMPFIND (see Fig.Q] 
middle panel; Fig. [2}. Nevertheless, the total flux of 728 Jy 
within these clumps, is nearly the same as for the CLUMPFIND 
clumps. 

2.2.3. SExtractor 

The third algorithm w e used to analyse th e data i s the SExtractor 
(SE) software from iBertin & Arnoutsl dl 996b . which origi- 
nally was developed to detect stars and galaxies in large 
surveys. However the algorithm has also been used to ex- 
tract sour ces in sub-mm obs e rvations of molec ular cloud 
structures (ISchuller et al.l 12009b ICoppin et all 12000b . Like the 
CLUMPFIND algorithm it uses the thresholding method to lo- 
cate the clumps. To deblend sources the algorithm contours them 
at 30 exponentially spaced levels between the peak values of 
the detected sources and the detection threshold and follows the 
structure downwards. At each of this levels it tests if there is an- 
other peak, with an intensity exceeding a certain fraction of the 
total intensity. If this is the case, it extracts a new source. For 
each pixel between maxima, SExtractor calculates its contribu- 
tion to each object by assuming a Gaussian profile, and converts 
this into a probability for that pixel to be assigned to a certain 
object. 
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Fig. 3. The relative difference of the GAUSSCLUMPS (green) 
and SExtractor (yellow) fluxes with respect to the CLUMPFIND 
fluxes of the 38 commonly identified clumps. 



SExtractor determined an rms noise level of about 
30 mJy/beam so we used a detection threshold of 2<x to make 
our results comparable to the other algorithms. Corresponding to 
our resolution we defined the minimal area of a detected clump 
to 10 pixel. SExtractor finds 432 sources (see Fig. Q] right panel; 
Fig. |3 with a total flux of 823 Jy, 72% of the total flux in the 
map. This is slightly higher than the results of the other algo- 
rithms. 

All the extracted sources of the SExtractor sample, as well 
as the clumps identified by the former algorithms, have fluxes 
above 1.6 Jy/beam, except 18 of them that we considered arte- 
facts and do not further take into account in our analysis. 

2.3. Common Sources 

The three used extraction algorithms are working along differ- 
ent principles. While CLUMPFIND assumes no specific clump 
profile, the other two algorithms presume a Gaussian shape. 
Furthermore the CLUMPFIND and SExtractor algorithms use 
thresholding methods, while the GAUSSCLUMPS algorithm 
works on residual maps. 

In order to check how well the extraction results of the differ- 
ent algorithms agree, we defined a sample of "common sources", 
where the same clump is detected by all three algorithms. In the 
crowded and filamentary clouds in the center of the Carina neb- 
ula, many clumps are blended; therefore, the different algorithms 
produce often quite different decomposition results. 

However, in the less crowded parts of our map, we could 
identify 38 rather isolated clumps which are detected by all three 
algorithms in a consistent way. For these, we compare the de- 
rived total fluxes in Fig. 3, where we plot the relative difference 
of the GAUSSCLUMPS and SExtractor fluxes with respect to 
the CLUMPFIND fluxes. The plot shows a generally reasonable 
agreement, although differences by factors up to ~ 2 - 3 are seen 
in the total fluxes for some of the faintest clumps in the com- 
mon sample. For the majority of the clumps the deviations are 
less than 50%, and for most of the brighter clumps very good 
agreement, with deviations of less than 10%, is found. From this 
we conclude that for sufficiently bright isolated clumps the three 
extraction algorithms yield consistent total fluxes. 



3. Results and Discussion 

3.1. Determination of clump temperatures and masses 

The masses of clumps can be calculated from the observed sub- 
mm flux of thermal dust emission, which is optically thin, via 



M = 



d F\ tot R 

B V (T d ) Ky 



(1) 



where F v>tot is the observed spectral flux integrated over the 
source and d is the distance of the molecular cloud. For the 
gas-to-dust mass ratio R and the d ust emissivity k v we used the 
values from ISchuller et all d2009) i.e. R = 100 and K^o^m = 
1.85 cm 2 g _1 . B v (T,j) is the Planck function for a given dust tem- 
perature T d (see also lPreibisch et af]|201 Id) . 

A difference in temperature of only a few degrees in the 
relevant regime can change the derived masses of the clumps 
by a factor of 2 - 3 and therefore strongly affect the derived 
C1MF of the region dStamatellos et al.ll2007l) . Hence, for a reli- 
able mass estimation a good approximation of the temperatures 
of the clumps is indispensable. 

Unfortunately, for single-wavelength data sets (such as our 
LABOCA map) no direct determination of the cloud temper- 
ature is possible. In many investigations it is assumed that all 
clou ds in a map would have a common temp erature, e.g. 15 K 
(e.g. iMookerjea et alj|200l iKirk et al]|2006l) . While this may 
be a reasonable approximation in the case of more or less iso- 
lated clouds in a quiescent environment, we believe that such 
an assumption is not valid in the case of the Carina Nebula. 
The clumps in different parts of the complex will be affected by 
strongly different levels of irradiation, depending on the physical 
distance from the nearest massive stars, and this should lead to 
considerable differences of the cloud temperatures. 

In general, clouds with high column density should be cooler 
than low-column density clouds, because their interior is bet- 
ter shielded from the external radiation field. Such a relation 
between column density and cloud temperature has been ob- 
served in several cases, and can provide us with importan t infor- 
mation for an estimation of cloud masses. iPeretto et al.l (1201 Ol) 
recently determined the temperatures and the column densities 
of a sample of 22 infrared dark clouds (IRDCs) from the Hi- 
GAL Galactic Plane survey (see iMolinari et al1l2010l) by ana- 
lyzing their spectral energy distributions. Their Fig. 6 shows the 
expected anti-correlation between column density and temper- 
ature. To quantify this relation, we used their column density 
pe ak data, performed an Ordinary Least Squares (OLS) Bisector 
fit dlsobe et al.l fl990) and found the relation 



log A^JcnT 2 ] = 25.6 - 0.22 log T d [K]. 



(2) 



Since the column density of a clump with a dust temperature 
T d is given by 



R 



B v (T d ) Q k v fim u 



(3) 



where F y max is the spectral peak flux, £2 is the beam solid an- 
gle, and ju is the mean molecular weight, we have two equations 
relating the observed flux, column density, and temperature. We 
can now solve for the temperature for which the column density 
computed from Eq. © is closest to the relation in Eq. (f2). 

This gives us two solutions, from which we can reject 
the first one as it leads to temperatures below 4 K. The sec- 
ond solution gives us, for all samples, temperatures between 
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Fig. 5. The temperatures and fluxes of the CLUMPFIND sample, overplotted on the LABOCA map. The color of the circles cor- 
responds to the calculated temperatures and the size of the circles to the measured flux (~ 2 — 373 Jy/beam). The position of 
r\ Car is marked by the black star. The black box represents the region shown in Fig. [1] centred around Tr 14 (RA(J2000)=161°, 
Dec (J2000)=-59°.6). 



8.5 - 18 .5 K, which are typical values for molecular clouds and 
clumps dBergin & Tafallall2007l) . Most of the clumps (CF: 92 %, 
GC: 86 %, SE: 87 %) have temperatures above 14 K (see Fig.gJ 
top panel). The column density distributions we derive for the 
different samples are shown in the bottom panel of Fig. |4] 

In Fig. [5] we show the location and the temperature of the 
clumps extracted with CLUMPFIND. The clump temperatures 
are shown by the color while the size of the circles increases 
with the measured flux F v tot (see also Fig. |2j». The clumps are 
assembled along filamentary structures and pillars. These pil- 
lars show cold dense clumps on their tips, within which star for- 
mation occurs. Most of the brightest clumps are located within 
the structures a few arcminutes to the west of the Tr 14 cluster, 
where the clouds are dense and massive. The ionizing radiation 
of the nearby cluster can only affect the cloud surfaces and does 
not pervade to the cool center of the cloud. 

3.2. The Clump Mass Function (CIMF) 

With these individual clump temperature estimates we now can 
derive the masses of individual clumps. We find masses between 
about 5M and 4.7 x 10 3 M o and the total mass extracted in 
clumps is about 42 x 10 3 M Q (see TableQ])- 



For the differential Clump Mass Function (CIMF) we plot 
the number of clumps dN within a logarithmic mass interval 
d\og(M) against the mass. Up to ~ 10 M Q our sample is not 
complete due to the mass sensitivity limit and the wide range of 
background levels at different locations in our map. For clump 
masses above ~ 10 — 20 M Q one can see the expected decrease 
in the number of clumps for rising clump mass, which roughly 
follows a power-law 



To derive the slope of the power-law tail of t he ma ss spectra, 
we use the method of iMaschberger & Kroupal (2009) to fit the 
distribution. They use a bias-free Maximum Likelihood (ML) 
estimator to determine the power- law exponent without binning. 
This method takes into account data only above a given lower 
cut-off limit to calculate the power-law index. In order to make 
sure that the fit is not affected by incompleteness effects at lower 
masses, we use a conservative lower cut-off limit at 50 M Q (gray 
dashed line in Fig. [6} for the fit. For the power-law index we find 
a CF = 1.95 ± 0.07, a GC = 1 .95 ± 0.06 and a SE = 1 -89 + 0.06 for 
the CLUMPFIND, GAUSSCLUMPS and SExtractor samples, 
respectively (see Table[TJ. 
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Fig. 6. The mass histograms derived with the three dif- 
ferent clump-finding algorithms CLUMPFIND (top), 
GAUSSCLUMPS (middle) and SExtractor (bottom). The 
black lines show the samples and their power-law slopes with 
the temperatures derived from the relation we got from Peretto 
et al. (2010). The dashed gray line shows the used lower cut-off 
limit for these samples. For the red line we assumed a constant 
temperature of 10 K for all clumps. If we assume a constant 
temperature of 20 K the slope does not change, but the mass 
distribution is shifted to lower masses (blue arrow). 



Fig. 4. The distribution of individual clump temperatures (top) 
and column densities (bottom), with statistical errors, as found 
from Eq. (0 and Eq. (|2]i for the three different clump-finding 
algorithms. 
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Table 1. The masses and power-law indices for the 
CLUMPFIND, GAUSSCLUMPS and SExtractor samples 



F M m M max M min M peak a* 
[10 3 M o ] [M Q ] [M Q ] [M e ] 

individual clump temperatures; lower cut-off limit: 50M o 
CF 687 422 4 652.2 52 15T8 1.95 ± 0.07 
GC 371 42.6 4 458.2 8.2 25.1 1.95 ± 0.06 
SE 414 48.2 2 972.8 5.7 15.8 1.89 ±0.06 

T=10K; lower cut-off limit: 1OOM 

CF 69.1 3 366.2 14.7 39.8 2.15 ±0.08 

GC 69.9 3 225.9 14.6 63.1 2.08 ± 0.06 

SE 78.8 2 676.9 15.4 39.8 1.99 ±0.06 

T=20K; lower cut-off limit: 32M a 

CF 2L0 1 024.4 43 123 2.14 ±0.08 

GC 21.3 981.6 4.5 20.0 2.06 ± 0.06 

SE 24.0 814.6 4.7 12.6 2.02 ± 0.06 



Notes. The results are shown for the individual temperature case and 
the two constant temperatures. For each temperature we use an accord- 
ing conservative lower cut-off limit. 
(fl) Number of clumps ( *> power-law index; dN/dM oc M~ a 




1 10 100 1000 10000 

M [M e ] 

Fig. 7. The mass distribution of the CLUMPFIND sample 
with its power-law slope (a = 1.95) in black. In green the 
GAUSSCLUMPS sample is over plotted with its log-normal fit. 
The dashed gray line shows the lower cut-off limit. 



3.3. CIMF assuming constant cloud temperature 

A difference in temperature of only a few degrees can change the 
derived masses of the clumps significantly. Here we investigate 
how the CIMF depends on the estimated temperature. For com- 
parison we therefore, also computed clump masses based on the 
assumption of isothermal clouds and considered two different 
values for the cloud temperature, 10 K and 20 K. These values 
are often assumed in sub-mm studies of molecular clouds ( e . g . 
Uohnstone et aT1l2000t iKirk et al.ll2006t ISchulleret al.ll2009h . If 
we assume constant temperatures for all clumps the overall 
shape of the CIMFs keeps the same, but we derive a steeper 
power-law slope (e.g. for 10 K acF = 2.15 ± 0.08, ace = 
2.08 + 0.06 and a SE = 1.99 + 0.06 for a lower cut-off limit of 
100M o ). The slopes for the two different constant temperatures 
stay the same within their errors, as a change of temperature cor- 
responds to a constant shift in mass (Fig. |6). The masses range 
between about 4 M G and 1 x 10 3 M Q for 20 K and about 14M 
and 3.3 x 10 3 M for a temperature of 10 K. The results are sum- 
marized in Table [TJ 

3.4. Power-law versus Log-normal 

iPeretto & FuTieil (120101) found, for a sample of gravitationally 
bound IPvDCs, a power-law index a ~ 1.8, while the slope 
of the mass spectra of unbound fragments steepens at the 
high mass end. They also show that the mass function of the 
fragments can be well described by a log-normal distribution. 
These results are in agreeme nt with the theoretical work of 
iHennebelle & Chabrierf (120081) and with studies of the probabil- 
ity distribution function of the column density (N-PDF) within 
molecular clouds which displays log-normal shapes for turbu- 
lent structures and power-law behaviour for gravitationally dom- 
inated clouds dKainulainen et alj201 UlBallesteros-Paredes et al.l 
1201 11) . 

When we calculate the power-law slope for different lower 
cut-off limits, starting from the peak mass to higher masses, we 
find that the power-law index for our CLUMPFIND sample is 
quite robust, within the errors, regardless of whether we assume 
individual or constant temperatures for the clumps. The slopes 
of the other two samples however become continuously steeper 
while shifting the cut-off point to higher masses. This indicates 



that GAUSSCLUMPS and SExtractor samples can be also de- 
scribed by a log-normal function as well as by a power-law 
while the CLUMPFIND sample shows a clear power-law behav- 
ior (Fig.|7J. We thus find that, the results of the before mentioned 
studies have to be used with caution while drawing conclusions 
on the physical state of a single molecular cloud, as the shape 
of the CIMF can highly depend on the clump extracting method. 
Whether the clumps are gravitationally bound or dominated by 
turbulence has to be confirmed otherwise, e.g. by determining 
their virial parameter dBertoldi & McKeel[l992b . 



4. Summary 

We analysed a large-scale sub-mm map of the Carina Nebula 
Complex to investigate the mole cular cloud structures of 
the region dPreibisch et al.l 1201 Id) . From this map we ex- 
tract the clumps with three common clump-finding algorithms 
(CLUMPFIND, GAUSSCLUMPS and SExtractor) to derive the 
CIMF. We therefore are able to compare the resulting CIMFs and 
to test the influence of the extraction method and the effects of 
the assumed temperatures on the CIMF (see Table [1}. 

The assumption of a uniform temperature for the whole area 
disregards that the different regions within the Carina Nebula 
Complex are affected by highly variable levels of irradiation 
and stellar feedback. Hence we used an empirical calibration of 
the relation between column density and t emperature, which we 
adopted from data of IPeretto et al.l (1201 0l) . to estimate the tem- 
peratures and masses of the individual clumps. The clump tem- 
peratures for all extraction methods lie between 8.5 - 18.5 K, 
which are characteristic temperatures for molecular clumps. 

From the subsequent resulting masses we derived the CIMFs 
of this region. For masses above 50 M Q all samples can be 
described by a power-law with index a around 1.9, as de- 
fined in Eq. This is in good agreement with the re- 
sults of ot her studies of clump mass distributions in molecu- 
lar clouds dElmegreen & Falgaronelll996l iKramer et al. 1 H9981 
ISchneider & Brooksl 2004 ). The slope of the power- law we 
foun d is also quite simila r to the slopes for cluster mass functions 
(e.g. lLada & Ladall200l . If we assume a uniform temperature 
for all clumps (10 K / 20 K) the power-law slope gets steeper. 
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The GAUSSCLUMPS and SExtractor samples show a lack 
of clumps in the low-mass range, while the CLUMPFIND al- 
gorithm finds for this mass regime at least twice the num- 
ber of clumps. So even when the high mass end of the thus 
found CIMFs of the GAUSSCLUMPS and SExtractor samples 
show a power-law behaviour their overall shapes are better de- 
scribed by a log-normal function. The CLUMPFIND sample, 
however, clearly follows a power-law. Hence we find the shape 
of the C1MF highly dependent on the method used to extract the 
clumps. 

Theoretical models predict that unbound turbulent clouds 
should have a log-normal C1MF, whereas clouds contain- 
ing gravitationally bound clumps develop a power-law tail 
dKainulainen et alJl201 lb iPeretto & Fulled l2010l) . In this sense, 
the shape of the C1MF contains very important information 
about the evolutionary state of a molecular cloud. However, 
our results show that a reliable determination of the shape of 
the C1MF (log-normal or power-law) is often difficult as it can 
depend on the specific source extraction algorithm. Therefore, 
conclusions drawn from a C1MF derived from a single extrac- 
tion algorithm should be taken with care to avoid an over- 
interpretation. 
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